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Abstract 

We study the combinatorial complexity of D-dimensional polyhedra denned as the 
intersection of n halfspaces, with the property that the highest dimension of any 
bounded face is much smaller than D. We show that, if d is the maximum dimension 
of a bounded face, then the number of vertices of the polyhedron is 0(n d ) and the total 
number of bounded faces of the polyhedron is 0(n d ). For inputs in general position 
the number of bounded faces is 0(n d ). For any fixed d, we show how to compute the 
set of all vertices, how to determine the maximum dimension of a bounded face of 
the polyhedron, and how to compute the set of bounded faces in polynomial time, by 
solving a polynomial number of linear programs. 



1 Introduction 



Bounds on the complexity of halfspace intersections and convex hulls 24 
and algorithms for constructing halfspace intersections and convex hulls 
[6j[7j[l4] 1 |22 , 27, 30|, have long been a mainstay in discrete and computational 
geometry. However, because of their inherent exponential dependence on the 
dimension of the input, these worst case bounds are only useful for polytopes 
of low to moderate dimension. In contrast, linear programming allows single 
vertices of halfspace intersections to be found numerically in polynomial time 



even for high-dimensional inputs 20 21 29 . Is there an intermediate family of 



polytope construction problems, with low complexity even for high-dimensional 
input, but capable of describing more complex sets of features than a single 
polytope vertex? 

Given the obstacles that have been identified to finding a truly output- 
sensitive halfspace intersection algorithm [l], it may not be sufficient to ask for 
an output-sensitive algorithm that works efficiently whenever an input polytope 
has low combinatorial complexity. Instead, in this paper, we identify a different 
parameter of instance complexity that does lead both to low complexity and effi- 
cient algorithms: the dimension of the bounded faces of a halfspace intersection. 
Figure [l] shows an example polyhedron whose bounded faces have been marked. 
These faces form a combinatorial structure that can be described equivalently 
in several ways: 
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Fig. 1: An unbounded 3-dimensional polyhedron, with a bounded subcomplex consist- 
ing of four 2-dimensional and three 1-dimensional maximal faces. On the right 
we show the projection onto a plane that is perpendicular to the unbounded 
direction. 

• The bounded subcomplex of an intersection of n halfspaces is the cell com- 
plex formed by the set of bounded faces of the intersection. 

• The subset of faces of a convex polytope that lie above some linear thresh- 
old is equivalent to a bounded subcomplex under a projective transforma- 
tion that takes the threshold hyperplane to the hyperplane at infinity. A 
transformation of this type necessarily transforms the faces of the polytope 
that are entirely above the threshold into bounded faces of the transformed 
polyhedron, but the faces that cross the threshold (if any exist) become 
unbounded. 

• The subset of faces of a convex polytope that are disjoint from some spec- 
ified facet is again equivalent to a thresholded intersection or a bounded 
subcomplex. To see this, one may consider a linear threshold determined 
by a linear function that is zero on the specified facet, positive on the 
points of the polytope disjoint from the facet, and negative on the other 
side of the facet. 

• The complexes considered above are equivalent under projective duality 
to the subset of faces of a convex polytope that are nonadjacent to some 
specified vertex. 

We show that, when the bounded subcomplex consists only of faces of dimension 
at most d, where d is significantly smaller than the ambient dimension, then 
it has smaller combinatorial complexity than the bounds given by the upper 
bound theorem. Specifically, it has 0(n d ) vertices, and 0(n d ) faces. More 
strongly, when the halfspaces forming the intersection are in general position, 
we obtain a tighter bound of 0(n d ) on the number of faces as well. Based on 
these combinatorial bounds, we provide algorithms for listing all the vertices 
or faces of the bounded subcomplex in polynomial time whenever d is a fixed 
constant. 
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2 Motivating Applications and Related Work 

One motivation for studying bounded subcomplexes comes from the tight span 
construction, a canonical method of embedding any metric space into a contin- 
uous space with properties similar to those of L°° spaces [8 13 19 . One way 



of defining the tight span, for a finite metric space with n points Pi, is to coor- 
dinatizc n-dimcnsional L°° space by n variables Xi, and to define a polyhedral 
subset of the space by the Q) linear inequalities 

Xi + xj > dist(pi,pj) 

for each possible pair Then, the tight span is the bounded subcomplex 

of this polyhedron. For metric spaces satisfying an appropriate general position 
assumption, the dimension of the bounded subcomplex is between |~n/3] and 
[n/2\ 1 10] , but certain combinatorially defined metrics, such as the metrics of 
distances on certain classes of planar graphs, can have tight spans of much lower 
dimension [15| . Our results bound the complexity of these low-dimensional tight 
spans and allow them to be constructed efficiently, generalizing our previous 
algorithms for constructing tight spans when they are homeomorphic to subsets 



of the plane 16 



Our results can also be interpreted as statements about the complexity of 
Delaunay triangulations for inputs satisfying strong convex position assump- 
tions. Delaunay triangulations are closely related to convex hulls: the Delaunay 
triangulation is combinatorially equivalent to the convex hull of a point set lifted 
to a sphere in one higher dimension, augmented by an extra point at the pole 
of the sphere, followed by the removal of all faces incident to that pole [4]. If a 
D-dimensional set of n points has the property that every interior point of the 
convex hull of the set belongs to a Delaunay triangulation feature of dimension 
at least D — d, then our results imply via this lifting relation that the Delaunay 
triangulation has 0(n d ) D-dimensional simplices. For instance, if the convex 
hull of a point set is a stacked polytope (a convex figure formed by gluing sim- 
plices facet-to-facet) and the Delaunay triangulation coincides with the gluing 
pattern of the simplices, then every interior point of the convex hull belongs 
either to one of the simplices or to one of the glued facets, so d — 1; in this case, 
the number of simplices in the Delaunay triangulation is exactly n — D. 

Additionally, many combinatorial optimization problems such as shortest 
path trees, minimum spanning trees, bipartite minimum weight perfect match- 
ings, and network flows can be expressed as linear programs, and the vertices of 
the corresponding polyhedra that have the k smallest values (according to the 
linear objective function used for the problem) correspond to the best k solu- 
tions of these optimization problems. From this point of view, it is of interest 
to understand the combinatorial behavior of the subsets of polyhedron vertices 
and faces satisfying a linear threshold constraint, or equivalently of the bounded 
faces of polyhedra. As a somewhat trivial example (easy to prove directly with- 
out resorting to our polyhedral point of view), consider the bin packing problem 
in which the input consists of a set of items with weights u>i and a capacity W, 
and a solution consists of a subset of the items with total weight at most W. If 
Xi is an indicator variable for the inclusion of item i, then the solutions are the 
vertices of the hypercube < Xi < 1 that satisfy the linear constraint 

^2 W * X i < w - 
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The number of halfspaces defining the hypercube is 2n, and the dimension d of 
the bounded subcomplex of the hypercube determined by the linear constraint 
is the maximum number of items that can be packed into a single solution. 
Therefore, our results imply that the total number of solutions is 0(n d ). 

Bounded subcomplexes have also been investigated in other contexts. In 
tropical geometry, the tropical complex defined by a set of points is also equiv- 
alent to the set of bounded faces of a polytope defined from the points 11 . 
Develin [9] studied the bounded faces of a certain polytope arising from a prob- 
lem in algebra, and showed that they are all isomorphic to subpolytopes of 



permutohedra. Queyranne 26 used polytopes to model scheduling problems; 
the polytope defined by Queyranne has a unique bounded facet. In connection 
with the tight span application, Hirai [18] showed that the bounded faces of a 
halfspacc intersection form a contractable complex. Herrmann et al. 17 con- 
sider the problem of computing the bounded subcomplex of a polyhedron, but 
they do not bound the complexity of the complex. Their algorithms assume 
that the vertices and facets of the polyhedron are both already known, and are 
output-sensitive given this information. 



3 Definitions 



Define a polytope to be the convex hull CH(S) of a finite set S of points, and 
a polyhedron to be the intersection C\T of a finite set T of closed halfspaces. 
Every polytope is a polyhedron; a polyhedron is a polytope if and only if it is 
bounded. 

The faces of a polytope or polyhedron P are the sets of the form PDH where 
H is a closed halfspace whose boundary is disjoint from the relative interior of 
P. Every face of a polytope is also a polytope, and every face of a polyhedron 
is also a polyhedron. The dimension dim / of a face / is the dimension of its 
affine hull aff /. With these definitions, P itself is a face, as is the empty set. By 
convention the dimension of the empty set is —1. Two faces are incident if one 
is a subset of the other. A face / of a polytope CH(S) may be identified with 
the set / n S, and a face / of a polyhedron C\T may be identified with the set 
of halfspaces of T whose boundary hyperplanes contain /. This identification 
is one-to-one: each face is identified with a unique set. The subset relation 
between faces is either the same as or the reverse of the subset relation between 
the sets the faces are identified with, for polytopes and polyhedra respectively. 

A vertex of a polytope or polyhedron is a O-dimensional face; it must be a 
singleton set, containing a single point. An edge of a polytope or polyhedron is a 
bounded 1-dimensional face, consisting of a line segment connecting two vertices. 
Polyhedra may also have rays and lines, unbounded 1-dimensional faces with 
one or zero vertices on them respectively. If a polytope or polyhedron P has 
dimension D, then a facet of P is a (D — l)-dimensional face. If D is also the 
dimension of the space containing the points or halfspaces from which P was 
defined, then P is said to be full- dimensional. More generally we define a d-face 
of P to be a face of dimension d, so for instance a 1-face may be an edge, a ray, 
or a line. 

If v is a vertex of a polytope or polyhedron P, let ft, be a hyperplane that 
separates v from all other vertices and all unbounded regions of P, and define 
the link of v to be the intersection PDh. Figure [2] shows an example. Different 
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Fig. 2: The link of a vertex of P is a polytope of dimension one lower than P. 

choices of h will give geometrically different intersections, but they are all related 
to each other by affine transformations, so as a combinatorial polytope the link 
is well defined; however, in our proofs we will need to refer to the link as a 
geometric object, defined as above, rather than as a combinatorial object, so in 
each such usage we will need to verify that the choice of h does not affect the 
result. The faces of the link correspond one-for-one, in an incidence-preserving 
way, with the faces of P that are incident to v. 

A full-dimensional polytope or polyhedron is simple if every vertex is incident 
to at most D edges. Equivalently, a polytope or polyhedron is simple if the link 
of every vertex is a simplex. A polyhedron defined as the intersection of a set 
of halfspaces in general position (meaning that a small perturbation of any of 
the halfspaces does not change the combinatorial structure of the polyhedron) 
is necessarily a simple polyhedron. Dually, a polytope is simplicial if all of its 
faces except for the polytope itself are simplices. A polytope defined as the 
convex hull of a set of points in general position is necessarily simplicial. 

If P is a polytope containing the origin, the polar polytope of P is the inter- 
section of a system of halfspaces, one for each vertex of P: if v is a vertex of 
P then the corresponding halfplane is the set of points w such that v • w < 1. 
Although defined as an intersection of halfspaces, the polar polytope is bounded 
and is therefore a polytope. The faces of the polar polytope correspond one-for- 
one with the faces of P but in a dimension-preserving way: an i-dimensional 
face of P corresponds to a (D — i — l)-dimensional face of the polar polytope. 
The polar polytope of the polar polytope is P again. If P is simple, its polar 
polytope is simplicial, and if P is simplicial, its polar polytope is simple. 

If P is a polyhedron or polytope and £ is a linear function, we define max^ P 
and min^ P to be the face on which I takes its minimum or maximum. If the 
minimum or maximum is unbounded, then we define max^ P or min^ P to be 
the empty set. A polyhedron is pointed if it has at least one vertex. If £° is a 
tangent hyperplane to a polyhedron P at a given vertex v, that is not tangent 
to any higher-dimensional face, and I is a linear function that is positive on the 
side of £° containing P and negative on the other side, then v = min^ P, so a 
polyhedron P is pointed if and only if there exists a linear function I for which 
min£ P is a single vertex. 

A polyhedral complex is a finite set C of polyhedra, all in the same ambient 
space, such that C contains each face of each polyhedron in C and such that the 
intersection of any two polyhedra in C is a face of both. We define a polytopal 
complex to be a polyhedral complex in which each member of C is a polytope. 

If P is a polytope or polyhedron, £ is a linear function, and B is any real 
number, let P e<B be the polyhedral complex formed by the set of faces of P 
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Fig. 3: Trees that may be realized as the bounded faces of two- and three- 
dimensional polyedra. 

such that, for every point p of a face in P i<B , £(p) < B. Analogously, for 
B = oo, define P e< °° to be the set of faces of P on which I is bounded. (This 
notation directly implies only that I is bounded from above, but we define it 
to mean that it is bounded from below as well.) We define dimP i<B to be the 
maximum dimension of a face in P i<B . 



4 Examples 

We begin with some examples of polyhedra with low-dimensional bounded faces. 
These are not so much intended as motivation for the bounded-face problem but 
rather as illustrations of the sort of behavior these polyhedra may have; for more 
in-depth motivating examples, see Section [2] 

• Let P be any (D — l)-dimensional polytope, having facets determined by 
the linear inequalities on ■ x > hi. Then the D-dimensional polyhedron 
defined by the linear inequalities • x — biy > is an unbounded cone 
over P: it has P as its cross-section (say in the hyperplane y = 1) but has 
only a single vertex at the origin. Examples of this form can be used to 
show that, even when dimP^ <B is small, the number of unbounded faces 
can be large. For instance, if P is a hypercube, defined as the intersection 
of n = 2(D — 1) facets, then the cone over P has 3™/ 2 + 2 faces despite 
having only a single bounded vertex. 

• For d = 1 and D = 2, an unbounded two-dimensional polyhedron with n 
facets always has exactly n — 1 vertices (Figure [3]). For d = 1 and D = 3, 
if T is a tree forming the boundary edges and rays of a partition of the 
plane into n unbounded polygons [5] , with exactly three polygons meeting 
at each vertex, then T may be lifted to an unbounded three-dimensional 
polyhedron with n facets and n — 2 vertices (Figure [3]). 

• More generally, a D-dimensional polyhedron may have any free tree T 
with vertex degree at most D as its complex of bounded faces. To see 
this, let P be a (D — l)-dimensional stacked polytope, the union of a set 
of simplices glued facet-to-facet with one simplex for each vertex of T and 
with two simplices glued together whenever the corresponding vertices of 
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Fig. 4: The bounded subcomplex of the Voronoi diagram of n points on the 
three-dimensional moment curve, for 4 < n < 7. 



T are adjacent. It can be shown by induction on the size of T that there is 
always a way of choosing simplices to glue in this way such that their union 
is convex. P can be lifted to a stacked polytope P + in D dimensions, by 
adding the same new vertex v to every simplex of P; again, P + can be 
constructed by gluing simplices in such a way that their union is convex. 
Let Q + be the polar polytope of P + , let Q be the polyhedron formed by 
intersecting all but one of the facet halfspaces of Q + (omitting the facet 
dual to v), and let i be a linear function that is zero on the facet of Q + 
dual to v and negative on Q. Then Q l<0 is isomorphic to T. Note that, 
in defining P, all points may be chosen to be in general position, from 
which it follows that in Q, all halfspaces may be chosen to be in general 
position. 

Any n points on the three-dimensional m ome nt curve (i,i 2 ,i 3 ) determine 
a Voronoi diagram with ("j 2 ) vertices [l2| . A standard lifting trans- 
formation takes this Voronoi diagram to an unbounded four-dimensional 
polyhedron, with n defining halfspaces. In this polyhedron, all facets are 
unbounded, so the maximum dimension of a bounded face is two. There 
are quadratically many bounded faces (Figure Eh. As this example shows, 
a nontrivial bound on d can still lead to examples with superlinear com- 
plexity. 

Let P be a four-dimensional pyramid over an octahedron, modified by 
folding its octahedral facet so that it is split into two square pyramids. 
Then the subset of faces of P that are nonadjacent to the apex of P are 
dual to a one-dimensional complex that takes the form of a single line 
segment: the two square pyramid facets are dual to the endpoints of the 
line segment, and the square ridge separating these two facets is dual to 
the line segment itself. However, if the vertices of P are perturbed into 
general position, then the square becomes a flat tetrahedral facet while 
each of the two square pyramids is split into two tetrahedra; the dual 
complex becomes two triangles joined at a vertex. This example shows 
that some instances cannot be perturbed into general position without 
increasing the dimension of the bounded subcomplex. 

Let Qd be the D-dimensional hypercube, formed by intersecting n = 
2D halfspaces determined by the inequalities < x% and x^ < 1, and 
let a be the linear function J2 x i- Then, for any integer < d < D, 
dim Q a p D ~ d ~ 1 / 2 = d. The number of vertices of Q < £f >D ~ d ~ 1 ^ 2 is J2i=o ( D ) ■ 
For any even integer n, and any constant value of d, this provides an 
example of a bounded subcomplex in general position with <d(n d ) vertices. 



5 Bounding the number of vertices 



8 




Fig. 5: No matter how we slice a tetrahedron with a plane, there is always a pair of 
faces (one of which could be the empty set) on both sides of the plane such 
that the affine hull of these faces is the full space. Lemma [T] claims that this 
is still true for any polytope, in any dimension. 



5 Bounding the number of vertices 

In order to help prove our main result, we begin with a technical lemma stating 
that if we slice a polytope by a hyperplane in such a way that one side of the slice 
contains only faces of low dimension, then the other side of the slice contains a 
face of high dimension. Figure [5] illustrates the lemma. 

Lemma 1: Let P be a polytope, and A be a hyperplane that does not pass 
through any vertex of P. Then there exist faces / + and /~ of P, on opposite 
sides of A, such that aff (/+ U /") = aff P. 

Proof: We use induction on dimension. Note that the empty set is also a 
face of P, so A need not intersect P. Let £ be a linear function that is zero on 
A, let P + = P i>0 be the complex of faces of P on which £ is entirely positive, 
and let P~ = P l<0 be the complex of faces of P on which I is entirely negative. 
As a base case, if P + is empty, then we may take / + = and /~ = P. As a 
second base case, if P + is a single vertex v, then we may take f + — v and f~ 
to be any facet of P disjoint from v; such a facet must exist, for otherwise P 
would be unbounded. 

We may assume without loss of generality (by perturbing A, if necessary) 
that no hyperplane parallel to A passes through two or more vertices of P. By 
the assumption of general position on A, £ does not take the same value on any 
two vertices of P. Let v be the vertex in P + minimizing £, and let L be the link 
of v (with an arbitrary choice of the intersecting hyperplane fi defining the link) . 
Let L + — L l>l ( v ) and L~ = L e< *( v h By the induction hypothesis, applied to 
the polytope L, we may find faces g + and g~ in L + and L~ respectively, such 
that aff (g + U g~) — aff L. (When P is full-dimensional, this affine hull is just 
the hyperplane /z.) The face g + of L + corresponds to a face /+ of P in which 
all points have a value of £ larger than or equal to £(v); since £{v) is strictly 
positive, / + belongs to P + . However, the face g~ of L~ corresponds to a face 
h of P such that v lies in P + and the remaining vertices of h lie in P~ . Let f~ 
be any facet of h disjoint from v, as in the second base case. 

Then 

ff + = (/+nM)c/+Caff(/+U/-). 
Additionally, v € / + , so 

g- = (hHfj,) C h C aff(uU/ - ) C aff(/+ U /"), 

and 

aff P = aff(LU v) = aff( 3 + U g~ U v) C aff(/+ U /"). 
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Fig. 6: A vertex uofa 3-dimensional polyhedron has at least three incident edges. On 
the left, two edges point "up" (in the direction of £), so v is the minimum of 
the face between them: IS] = 1. In this case, the only face that is necessarily 
bounded is the edge that points down (which has dimension 1). On the right, 
only one edge points up, so v is the minimum of the intersection of the two 
faces on both sides of it: |S| = 2. In this case, however, at least two edges point 
down, so the face between them is necessarily bounded (which has dimension 
2). 

But / + U / _ is a subset of P, so its affine hull cannot be a proper superset of 
the affine hull of P, and the two affine hulls are equal, as desired. M 
As we now show, the vertices of a polyhedron with low-dimensional bounded 
faces can be associated to small sets of facets of the polyhedron, leading to poly- 
nomial bounds on the number of vertices of the polyhedron. Figure [6] illustrates 
the lemma for 3-dimensional polyhedra. 

Lemma 2: Let P be a full-dimensional pointed polyhedron, defined as the inter- 
section of n D-dimensional halfspaces. Let £ be a linear function with a vertex of 
P as its minimum, and suppose that £ does not take the same value at any two 
different vertices of P. Let B be a given real number, and let d = dimP e<B . Let 
v be any vertex of P e<B . Then there exists a set S of d of the halfspaces defin- 
ing P, such that v = min^(P n A), where A is the intersection of the boundary 
hyperplanes of S. 

Proof: Define a linear function t v by £ v (x) = £(x) — £(v). Then £ v differs 
from £ by a simple translation; therefore, it has the same maxima and minima 
as £ on each face of P. Let V be the link of v, let V + be the complex of faces of 
V on which £ v is positive, and let V~ be the complex of faces of V on which £ v 
is negative; this definition of V + and V~ does not depend on the intersecting 
hyperplane used to define the link geometrically. 

By Lemma [l] we can find faces / + and /~ in V + and V~ respectively, 
such that aff(/ + U /~) = aff V. Since dimF = dimP — 1 = D — 1, we have 
dim /+ + dim /~ > D — 2. Then these two faces / + and / _ correspond to faces 
g + and g~ in P, such that v = min^ g + = max£ i?~, and such that dim g + + 
dime/ - > D. But we know that dim<?~ < d by assumption, so dimg + > D — d. 
Therefore, g + may be represented as the intersection of P with a set of the 
boundary hyperplanes of at most d of the defining halfspaces of P. If A is the 
intersection of these hyperplanes, then v — min^ g + = min^Pn A) as desired. 

Theorem 3: Let P be a full-dimensional pointed polyhedron, defined as the in- 
tersection of n P-dimensional halfspaces, let £ be a linear function which attains 
a minimum on P, let B be given and greater than the minimum value of £, and 
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suppose that dimP^ < 
most 



Proof: Lemma [2] gives a formula v = mmg(P n A) by which each vertex 
v of P may be obtained from a subset of at most d of the defining halfspaces 
of P. We may augment any set of halfspaces determining v in this way so 
that it has exactly d halfspaces, by making an arbitrary choice of a sufficient 
number of the other halfspaces whose boundaries are incident to v. There are 
Q) sets of exactly d halfspaces, each one determines at most one vertex, and all 
vertices can be determined in this way, but some vertices may be duplicated. 
In particular the vertex min^ P belongs to at least D facets, and is the vertex 
that is determined by each of the (^) different d-tuples of these D halfspaces, 
giving rise to the correction terms in the formula. M 

Corollary 4: Let P be a polyhedron, defined as the intersection of n D-dimensional 
halfspaces, let £ be a linear function which attains a minimum on P, let B be 
given, and suppose that dim.P e<B = d < D. Then the number of vertices of 
P l<B is 0(n d ). 

Proof: We may assume without loss of generality that P is full-dimensional, 
for otherwise we may restrict our attention to the affine hull of P and eliminate 
any halfspaces that contain the entire affine hull. 

If P is not pointed or B does not exceed the minimum value of £, then 
P e<B has zero vertices and the result is obvious. Otherwise, we may ignore 
the — + 1 correction terms in Theorem [3J obtaining a formula which grows 
asymptotically as 0(n d ) and is independent of D. M 

Because of the upper bound theorem, these bounds can only be tight when 
d < D/2. However, as our main interest is for bounded d and unbounded D, this 
is not a significant limitation. The formula of Theorem [3] is tight when d = 0, 
and the examples of polyhedra with arbitrary .D-ary trees as their bounded 
subcomplex show that it is also tight for d = 1. The example of the Voronoi 
diagram of points on the moment curve shows that the asymptotic bound of 
Corollary [2] is also tight when d = 2 and D — 4, and more strongly that the 
formula of Theorem[3]is tight to within a 1 — 0(l/n) factor in this case. The hy- 
percube example d shows that the asymptotic bound of Corollary [4] 

is tight for any constant d. 



= d < D. Then the number of vertices of P t<B is at 




6 Euler's formula 

Letting Q(P e<B ) denote the set of faces of P e<B , define the Euler characteristic 

x(P e<B )= £ (-i) dim/ . 

fe®(P t<B ) 

Recall that we include the empty set as a face, as well as P itself in the case that 
maxi P < B. As is well known, with this definition, x(P) = f° r am7 convex 



polytope P; Ziegler 31 gives a nice proof of this fact based on shelling. 

As we now show, the Euler characteristic of P e<B behaves similarly to the 
Euler characteristic of an entire polytope. This fact will lead to improved bounds 
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on the number of faces of P e<B , because it will imply that the ei-dimensional 
faces can be charged against other faces of lower dimension. Our proof can be 
viewed as a dual form of the shelling proof of Euler's formula for polytopes. 

Theorem 5: Let P be a polytope or polyhedron, let £ be a linear function, and let 
B be either oo or a real number. Additionally, suppose that p i<B is nonempty. 
Then x{P e<B ) = 0. 

Proof: We use induction, both on the dimension of P e<B and on the number 
of vertices of P e<B . As a base case, the result is clearly true if P e<B consists 
only of a single vertex: for, in that case, F contains two faces, the empty set of 
dimension —1 and that vertex, and these two faces have opposite signs in the 
sum defining x(P £<B )> cancelling each other out. 

Otherwise, P l<B has at least two vertices. We may assume without loss of 
generality (by perturbing £ if necessary) that no two vertices of P e<B have the 
same value of £. Let v be the vertex of P e<B with the maximum value of £, and 
let B' = £{v). Then, by induction, x{P t<B ') = 0. 

pKB differs f rom pi<B ky get of faces incident to v. Let L be the link 
of v in P; then the face structure of L l<B does not depend on the intersecting 
hyperplane used to define the geometry of L, and the faces incident to v in 
pKB are - m one _to-one correspondence with the faces of L e<B . For instance, 
the empty face of L l<B corresponds to the vertex v of P e<B , and each vertex 
u of L e<B corresponds to an edge uv of P e<B . In this correspondence, the 
dimension of a face in L t<B is one less than the dimension of the correspond- 
ing face of P e<B . Therefore, each face of L e<B makes a contribution to the 
Euler characteristic of L e<B with the opposite sign to the contribution of the 
corresponding face in P e<B . Therefore, x{P i<B ) = x{P t<B ') - x(L e<B ')- But 
L l<B has lower dimension than P l<B , and is non-empty (it includes at least one 
vertex corresponding to an edge that can be reached from v by a single step of 
the simplex method), so by induction on dimension, x{L t<B ) = 0. Therefore, 

x(P e<B ) = x(P e<B ') - x(L e<B ') -o-o^o. m 

As a simpler proof for the special case in which P is a polytope and / is a 
facet of P on which £ = B, let e be sufficiently small that there are no vertices 
of P for which £ is between B and B — e, and let A be the hyperplane £ = B — e. 
Then 

x(P i<B ) = x(P) + x(P n A) - x (f) = o + o-o = o. 
7 Bounding the number of faces 

Theorem 6: Let P be a polyhedron, defined as the intersection of n D-dimensional 
halfspaces, let £ be a linear function which attains a minimum on P, let B be 
given, and suppose that dimP e<B = d < D. Then the number of faces of P l<B 
is 0(n d2 ). 

Proof: By Theorem [5] the number of rf-dimensional faces is no larger than 
the total number of faces of all lower dimensions. Therefore, we need only count 
faces of dimension up to d — 1. But each such face is the intersection of P with 
the affine hull of a set of at most d vertices of P e<B . There are 0(n d ) vertices, 
so there are 0(n d ) sets of at most d vertices, and therefore 0(n d ) faces. M 

For instances in general position, a considerably sharper bound may be ob- 
tained. 
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Theorem 7: Let P be a polyhedron, denned as the intersection of n D-dimensional 
halfspaces, let i be a linear function which attains a minimum on P, let B be 
given, and suppose that dimP e<B = d < D. Additionally, suppose that the 
halfspaces defining P are in general position and that P has iV vertices. Then 
the number of i-dimensional faces of P e<B is at most ATm. 

Proof: We charge each i-dimensional face / of P e<B to the vertex max^ /. 
Let v be any vertex, and let L be its link. Then the faces charged to v correspond 
one-to-one with the faces of L e<e ( v \ By the general position assumption, L is 
a simplex, and therefore L l<1 ^ is also a simplex. Because the faces of P e<B 
correspond to faces with dimension one less in L e<e ( v \ dim L e<e ^ < d — 1. 
Thus, the number of i-dimensional faces of P e<B that are charged to v is at 
most the number of (i — l)-dimcnsional faces of a (d — l)-dimensional simplex, 
which is (f) . m 

Corollary 8: Let P be a polyhedron, defined as the intersection of n £>-dimensional 
halfspaces, let £ be a linear function which attains a minimum on P, let B be 
given, and suppose that dimP f<B = d < D. Additionally, suppose that the 
halfspaces defining P are in general position. Then the number of faces of P l<B 
is 0(n d ). 

Proof: By Theorems [3] and [7J adding the bounds on the the number of 
faces of P e<B over all dimensions up to d gives a total that is less than (™) 2 d < 

„d2 d _ nt„d 



0{n d ). 



8 Algorithms 

We now discuss algorithms for finding the bounded subcomplex of a polyhe- 
dron. There are many alternative solutions available, depending on whether 
the dimension d of the subcomplex is known or unknown to the algorithm, on 
whether only the vertices of the polyhedron need to be found or whether the 
whole bounded subcomplex is to be constructed, and on whether we use as a 
subroutine the algorithm for computing bounded subcomplexes of Herrmann et 



al. 17 or whether we directly enumerate the faces of the bounded subcomplex, 
using linear programming to test whether each face is bounded. 

Our algorithms will necessarily involve the solution of linear programs and 
linear feasability problems. Especially in the case of inputs that are not in 
general position, it is important that these linear programs are solved in an exact 
model of computation that allows solution vertices to be compared for equality: 
it is possible that a small perturbation of an input problem, as might occur 
due to round-off error in a non-exact numerical linear programming algorithm, 
could significantly increase the dimension of the bounded subcomplex. However, 
this sort of exact computation model is standard in computational geometry 
algorithms, and solutions are available when the ambient dimension D is of 
moderate size [23] . Even when D is large, strongly polynomial algorithms are 
known for some special cases of linear programming, such as the case with 



two variables per inequality that arises in the tight span construction 25 . We 
let L denote the time to solve a linear program of ambient dimension D and n 
constraints, or to determine whether such a program is infeasible or unbounded. 

When the input consists of a set of n halfspaces together with the dimension 
d of the bounded subcomplex, we have the following results: 
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Fig. 7: Lemma [9] states that if there exists a vertex of P that is the minimum on 
the intersection of d hyperplanes but no set of less than d hyperplanes, then 
there also exists a vertex of P that is the minimum on the intersection of d — 1 
hyperplanes but no set of less than d — 1 hyperplanes. In the figure, d = 2. 



• We can construct all vertices of the given polyhedron in time 0(n d L). This 
follows immediately from Lemma [2] which gives a formula v = min^(Pn A) 
allowing each vertex v of P to be obtained as the solution to a linear 
program in a subspace of dimension D — d. 

• If the input is in general position, we can construct the bounded subcom- 
plex in time 0(n d L + n 5d ). This method generates the vertices as above, 
exhaustively tests each vertex-facet pair to find all the vertex-facet inci- 
dences, and then uses the algorithm of Herrmann et al. 17 to construct 
the bounded subcomplex from the vertex-facet incidences. The algorithm 
of Herrmann et al. takes time 0(N 2 f3cf> + N 3 <fr 2 ), where N is the number 
of vertices of the given polyhedron (here at most 0(n d )), (3 is the number 
of vertex-facet incidences (at most 0(n d+1 ), and (j> is the number of faces 
of the bounded subcomplex (again, at most 0(n d )). 

• If the input is not in general position, the same method takes time 0(n d L+ 

n 2d 2 +3d^ 

• If the input is in general position, we can alternatively construct the 
bounded subcomplex in time 0(n 2d L). The method is to maintain a list 
of faces of the bounded subcomplex, initially containing all the vertices. 
Then, for each face / added to the list, and each vertex v, we compute 
the affine hull of / U v, use linear programming subproblems to determine 
whether this affine hull lies on the boundary of the polyhedron and has a 
vertex of maximum value, and if so check that the face contained in this 
affine hull is distinct from the ones already discovered. Whenever we find 
a new face in this way we add it to the list. There are 0(n 2d ) face-vertex 
pairs, each taking time O(L) to check, so the time bound is as given. 

• If the input is not in general position, the same method takes time 0(n d +d L). 

When d is not given as input (as seems more likely to occur in most appli- 
cations of this problem) , constructing the bounded subcomplex becomes more 
complicated. It is possible to interleave the construction of vertices (as the solu- 
tions of linear programming subproblems with decreasing dimension) with the 
construction of faces of increasing dimension, but how can we tell when we have 
found everything? An answer is provided by the following lemma, which we 
illustrate in Figure [7j 

Lemma 9: Let H be a given set of halfspaces, let P be the D-dimensional poly- 
hedron formed as the intersection of the halfspaces in H, let d < D be a number, 
and let I be a linear function. For any given set S of halfspaces in H, let Ps 
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denote the intersection of P with the hyperplanes bounding the halfspaces in S. 
Let V be the set of all vertices of P of the form min P$ for sets S of at most d 
halfspaces in H, and F be a set of bounded faces of P that includes the empty 
face. Suppose additionally that the following three conditions are all true: 

1. For each set S of at most d + 1 halfspaces in H, ming P$ belongs to V. 

2. Every face in F has dimension at most d. 

3. For every pair (f,v) of a face / in F and a vertex v in V, with A being 
the affine hull of / U v, one of the following three possibilities is true: 

(a) A contains an interior point of P and d < D — 1, 

(b) Some face f £ F has affine hull A, or 

(c) I is unbounded on A n P. 

Then F is the bounded subcomplex of the given polyhedron. 

Proof: We assume for a contradiction that there is a bounded face / that 
does not belong to P; among all such faces, let / have the minimum possible 
dimension, and let v = max^ /. If / had dimension at most d, then by Lemma[2] 
v would belong to V, and the pair (f',v) would violate condition (3) of the 
lemma where /' is any facet of / nonincident to v. We can assume without loss 
of generality that the lower link of v (the subset of the link consisting of the faces 
whose values in I are entirely less than that of the vertex itself) has dimension 
exactly d: it has dimension at least d, because it contains a face of dimension 
dim / — 1 corresponding to /, and if it had a higher dimensional face then we 
could reduce the dimension by one unit per step by moving from / to a facet of 
this higher dimensional face; condition (2) of the lemma ensures that, at each 
step, we continue to have a face / that does not belong to F . But, if the lower 
link of v has dimension exactly d, then by Lemma [2j v would be representable 
as min^ Pg for a set S of at most d + 1 halfspaces of H, violating condition (1) 
of the lemma. This contradiction shows that no such face / exists. M 

Theorem 10: Given a polyhedron whose bounded complex has dimension d (with 
d unknown to the algorithm) we can construct the bounded complex in time 
0(n d +d L). If the input is in general position, we can construct the bounded 
complex in time 0{n 2d L). 

Proof: We apply the algorithm for known d, for increasing values of d, until 
the condition of the lemma is met. Testing the condition takes time 0(n d +d L), 
or 0{n 2d L) for inputs in general position, matching the time for the known-c? 
algorithm on the correct value of d. M 

It would be of interest to determine whether the method of Herrmann et 



al. 17 , which uses many fewer linear programs at the expense of greater running 
time in the other parts of the algorithm, can be adapted to the case where d 
is unknown. In all cases the solution is found within an amount of time and 
a number of linear programming subproblems that is polynomial for any fixed 
value of d. 
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9 Discussion 

We have shown that, for any fixed d, bounded subcomplexes of polyhedra that 
have dimension at most d have polynomial complexity. Our bounds on the 
numbers of vertices of the bounded subcomplex are tight for small values of d. 
However, our bounds on the numbers of faces of higher dimensions do not ap- 
pear to be tight. In other bounds on the complexity of polytopes, it is possible 
to considerably simplify the problem by assuming that the input is in general 
position; for our problem, too, such an assumption would be very helpful (it 
would lower the exponent of the polynomial from d 2 to d), but it is not al- 
ways possible to perturb the input into general position without changing its 
dimension. Nevertheless, we would like to know whether instances that are not 
in general position have the same complexity as, or higher complexity than, 
instances that are in general position. 

For the problem of constructing tight spans, our results are even more unsat- 
isfactory, because the number of halfspaces is already quadratic in the number 
of points of the input metric space. For instance, when d = 2, we get a bound 
of 0(n 8 ) on the complexity of the tight span of an n-point metric space. This 
stands in contrast with the bounds from another paper of 0(n) on the combina- 
torial complexity of the tight span and 0(n 2 ) time to construct it (optimal since 
the input distance matrix has size 0(n 2 )) under a stronger two-dimensionality 



assumption, that the tight span is homeomorphic to a subset of the plane 16 
Obtaining tighter bounds for c?-dimensional tight spans, or even for the case 
d = 2 without the requirement that the tight span form a planar set, would be 
of interest. 
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